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Abstract 

We continue here the study of free extreme values begun in [3]. 
We study the convergence of the free point processes associated with 
free extreme values to a free Poisson random measure ( [15] ,[!])■ We 
relate this convergence to the free extremal laws introduced in [3] 
and give the limit laws for free order statistics. 

1 Introduction 

In classical probability theory, the theory of extreme values for i.i.d. random 
variables is elementary and well understood. Recently, a similar theory has 
been introduced in the context of free probability theory, in which the role 
of independent random variables is played by freely independent operators in a 
Hilbert space ([3]). The asymptotic behavior of the maximum of N free operators 
is given in [5], where the maximum is taken for the spectral order relation on 
operators ([1], [13]). The theory emerging is then parallel to the classical theory 
for maxima of i.i.d. random variables. In this paper, we make the next step in 
developing this parallel picture. We study the behavior of the full point process 
of normalized free extreme values. We show that it converges to a free Poisson 
random measure as soon as the normalized free maximum converges. One should 
notice that the notion of "free order statistics" is not readily available. Indeed, 
the notion of a "second largest statistic" is not at all clear. This difficulty is 
mirrored in the nature of the limiting object. Free Poisson random variables are 
not discrete. We will see (in Theorem [3]) that our main convergence theorem 
(Theorem [IJ leads to results with no classical analogs for order statistics. 
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The basic element in both classical and free theory of extremes is a probabil- 
ity measure /x. In the classical case, we take a sequence of i.i.d. random variables 
Xi, distributed according to /Lt, and introduce their order statistics, i.e., order 
them in non increasing order: 

X(0) >x(i) >X(2) > ... >x("-i), 

so that is the maximum of the n-sample, X'^^^ the second largest value 

and so on. The basic question is to describe the asymptotic behavior of the 
distribution of these order statistics once properly normalized, when n tends to 
oo. 

Let Fn^k denotes the distribution function of the normalized order statistics 
— - — -J for well chosen normalization constants a„ and 6„ 

< t 



Fn,k{t)=P 

The first question addresses the behavior of the maximum, i.e., the asymptotic 
behavior of F„.o. It was shown in the classical works by [8], [7], and [9] that 
there are only three types of possible limit laws, to which Fn^ can weakly 
converge. These laws (Weibull, Frechet or Gumbel) are called "extreme value 
distributions" : 



Type I : G (x) — exp (— e ^) , — oo < a; < oo 
Type II : G{x)-- 



0, x < 0, 

exp (— , for some a > 0, x > Q; 



Type III : G {x) 



exp (— (— a:)") , for some a > 0, a; < 0, 
1, a; > 0. 



Moreover, the nature of the max-domain of attraction of these extreme value 
distributions is well known as well as the possible choices for normalization 
constants ([IT], [H]). 

In the free probability context, a sequence of free self-adjoint operators Xi 
is taken, such that each of Xi has the spectral probability distribution ji. In 
recent work a maximum operation was defined which maps any n-tuple 
of self-adjoint operators to another self-adjoint operator, which is called their 
maximum. The definition is based on the so-called spectral order for self-adjoint 
operators: A < B iff all spectral projections l(_cx),t] (^) are greater than or 
equal to the corresponding spectral projections \(-oo,t] (B) ■ 

The spectral order is stronger than the usual order on operators, according 
to which A < _B iff i? — A is non-negative definite. The main benefit of the 
spectral order is that the set of all self-adjoint operators forms a lattice with 
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respect to this order. In particular, if S is the set of all operators C such that 
Ai ^ C for each of A i, ... , An, then S has a unique minimal element which is 
called max {^i, . . . , An} ■ This property does not hold if self-adjoint operators 
are considered with respect to the usual order on operators. Note, however, 
that the lattice of selfadjoint operators with respect to the spectral order is not 
a vector lattice in the sense that A — B ^ does not imply that A'^ B. For a 
counter-example and other information about the spectral order, see ^T3\ . 

By analogy with the classical case, the sequence of normalized maxima is 
defined as 

max {(Xj - bnl) /an} 

l<i<n 

where the maximum here is understood with respect to the spectral order. Then, 
^n'o'^ (x) is defined as the spectral distribution function of this normalized max- 
imum. 

In [3] the following question is solved: When does the sequence of F^J'^'^ 
converges weakly? 

The answer to this question is very similar to the answer in the classical 
case: There are only three possible types of limit laws, and for a given /i, the 
distributions fJJ'q'^ can converge to only one of them: 



Type I : G^'''"' (x) = | ^ 
II : Gf-""" (x) = I ^ 



0, a; < 0, 

e~^, a; > 0; 



Type 



X <l, 

for some a > 0, a: > 1; 



], a;<-l 
Type III : G^'''''' {x) ^ { 1 - (-x)" , for some a > 0, -1< x < 0, 

L, a;>0. 

As in the classical case, this allows defining domains of attraction of the free limit 
laws. Similar to the results about sums of free operators ([5]), an important fact 
is that, even though the limit laws are different in the classical and free cases, 
the domains of attraction are the same as well as the normalization constants! 
More precisely Fn^ converges weakly to the extreme value distribution G[x) iff 
-^n cT'^ converges weakly to G^^^'^ of the same type as G {x) . 

This rigid link between classical and free probability theory for extreme 
values is thus exactly similar to the analogous results for sums of i.i.d. random 
variables, as developed in ([6]). 

In order to investigate this situation further, let us return to the classical 
case and consider the random point process 



Nn=J2S^ 



{Xi-b„)/a„ 
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The next question of classical extreme value theory is to understand the conver- 
gence of this point process. This question is naturally related to the convergence 
of the distributions Fn_k- If M is in the domain of attraction of a classical ex- 
treme value distribution G{x), or equivalently if Fnfl converges to G{x) for 
some choice of normalization constants a„ and 6„, then the point process Nn 
weakly converges to a Poisson random measure with intensity measure A (dx) 
with X{x,oo) = —\ogG{x). Conversely, if 7V„ weakly converges to a Poisson 
random measure with the intensity measure A (dx) , then the distribution of any 
order statistics Fn^k converges to a limit law Gj^,) which is easily computable 
from X(dx) or equivalently from G(x), see below or ( [1 4J ) . 

What is the free analogue of the point process -/V„? To motivate our defi- 
nition, note that we can think about Nn as a linear functional on the space of 
bounded measurable functions: {Nn, /) —'■ X^iLi / ii-^i ~ Z*^") • This func- 
tional takes values in the space of bounded random variables. We will define a 
free point process analogously. We begin with a slightly greater generality and 
associate a free random process to any triangular array of free random variables. 

Let A be the set of densely- defined closed operators affiliated with a von 
Neumann algebra A, and let Boo (R) denote the set of all bounded, Borel mea- 
surable functions / : R ^ M. 

Definition 1 Let Xi.„ G -4, (i — 1, . . . ,n; n — I, . . .) be a triangular array 
of freely independent self-adjoint variables. Then the free point process Mn 
associated with the array Xi^n *s the sequence of A-valued functionals on Boo (R); 
defined by the following formula: 

n 

{MnJ) :=^/(X„0. 

i=l 

The triangular array of free variables that we use in applications to free 
extremes is, of course, Xi.„ = (Xi — &„) /a„, where Xi is a sequence of free 
self-adjoint variables. 

We can also define the concept of weak convergence of a free point process 
as a weak-* convergence of the corresponding functionals. In the classical case, 
after a suitable scaling, the point process Nn converges to a Poisson random 
measure. It turns out that in the non-commutative case the free point process 
converges to a free Poisson random measure, which was recently defined in [15j 
and The following three theorems are the main results of our paper. 

Theorem 1 Let G{x) be a classical extreme value distribution, i.e. a Gumbel, 
Frechet or Weibull distribution. Let x = inf {x : G (x) > 0} and define a mea- 
sure X{dx) on \x,oo) by the equality A((a;, oo)) — — \ogG{x). The following 
statements are equivalent: 

(i) fi belongs to the domain of attraction of the classical extremal limit law 
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G{x), i.e., for some constants an and bn the distribution Fnfl converges weakly 
to G{x); 

(a) fJ, belongs to the domain of attraction of the free extremal limit law G-^^'^'^ , 
i.e., for some constants a„ and bn the spectral distribution of the normalized 
free maximum, fJ^^q'^ converges weakly to G^'^^'^; 

(Hi) For some a„ and bn, the point process Nn weakly converges on (x, oo) to 
the Poisson random measure with intensity A {dx) ; 

(iv) For some On and bn, the free point process M„ weakly converges on [x, oo) 
to the free Poisson random measure with intensity A {dx) . 

In case one of the equivalent conditions in Theorem [1] is satisfied, then all 
the normalisation constants a„ and &„ can be taken to be the same in all four 
statements. 

The equivalences of (i) and (iii) follows from the results in [T3] (see, e.g.. 
Section 4.2.2 on page 209), and the equivalence of (i) and (ii) was proved in [5]. 
Thus, we only need to prove the equivalence of (i) and (iv). 

The equivalence of (i) and (iv) will be seen, in Section [31 as a consequence 
of the following more general result about convergence of free point processes. 
Recall that a measure is called Radon if /i [K] < oo for every compact K. 

Theorem 2 Let Xi^n be a triangular array of free, self-adjoint random variables 
and let the spectral probability measure of Xi n be fin- Let X be a Radon measure 
on D C R. The free point process M„ associated with the array Xi,n converges 
weakly on D to a free Poisson random measure M with the intensity measure A 
if and only if 

nfin (A) ^ A (A) (1) 

for every Borel set A <Z D. 

We now want to show what Theorem [1] implies for free order statistics. We 
begin by recalling basic facts about the classical theory of extreme values. If the 
measure n is in the domain of attraction of the extreme value distribution G{x), 
then as mentioned above, the convergence of the point process iV„ implies easily 
the convergence of order statistics. Indeed with the notations introduced above, 
it is easy to relate the distribution Fn.k of the normalized fc-th order statistics 
to the point process Nn, through the basic identity: 

= P[Nn{t,^) < k]=E[l[o^k]i{Nn,^t,oc)))] ■ 

This implies easily that the distribution Fn.k of the properly normalized 
order statistics weakly converges to the distribution 

k 

Gik){t) = 2^e . 



Fn,k (t) = P 



- bn 



< t 
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We now want to sec how this translates in the free context. More precisely, let 
Xi,...,Xn be freely independent self-adjoint variables with (possibly different) 
distribution functions Fj. Consider M„ the free point process associated with 
the sequence Xi and let 

n 

Yn {t) := (M„, = ^ {Xi) . 

1 



IS 



Definition 2 For every real k>0, we say that fJ^^^'^ {t) := E [l[o,/c] {Yn {t))] 
the distribution function of the k-th order statistic of the sequence Xi,. . . , X„, 
and that it is the k-th order free extremal convolution of the spectral distribution 
functions Fi. 

Note that the definition is valid not only for all integer k but also for all 
non- negative real k. 

One question that immediately arises is whether we can define an operator, 
for which the distribution F^'^^^ (t) would be a spectral distribution function? 
The answer to this question is positive. The condition t' > t implies that 
Yn {t') < Yn (t) and l[o,k] (^n {t')) > l[o,fe] {Yn {t)) . Therefore, as t grows, the 
operators l[o,fe] {Yn {t)) form an increasing family of projections and we can use 
this family to construct the required operator by the spectral resolution theorem. 

Definition 3 For every real k>0, let 

= J t dl[o,fe] {Yn {t)) . 

We call Z^^^ the k-th order statistic of the family X^. 

Prom the construction it is clear that Fn^k {t) is the spectral distribution 
function of the operator Z^^^ . 

In complete analogy with the classical case the limits of these free extremal 

convolutions can be computed using the limits of free point measures. If G {x) 
is one of the classical limit laws, then we use G^~^^ {x) to denote the functional 
inverse of G {x) . Let 



t- {k) 




to{k) 




t+ (k) 





exp 

1 

e 

exp 



Let A (i) 



■ logG (i) and pt (0 = (27rO"' ^4^ - (1 - A (i) + Cf- 
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Theorem 3 Suppose that measure fj, belongs to the domain of attraction of a 
(classical) limit law G {x) and a„, 6„ are the corresponding norming constants. 
Assume that Xi are free self-adjoint variables with the spectral probability mea- 
sure /i and let F^^^^ (t) denote the distribution of the k-th order statistic of the 
family {Xi — bn) /cLn, where i — Then, as n ^ oo, the distribution 

^rJk'^ (i) converges to a limit, [t) , which is given by the following formula: 



Fik) (t) 



0, ift<t^, 
1- A(t)l[o,i)(fc), ift>t+. 



It turns out that in the particular case of the 0-order free extremal convo- 
lutions, their limits coincide with the limits discovered in [3] (see Definition 6.8 
and Theorems 6.9 and 6.11): 

Fio)(t) = (1 - e-*) 1(0,00) W ; 

^(0) w = (i - ^) 1(1.-) w ; and 

F[^)'{t) = (l-|tr)l(_i.o)W + l[o,oo) W, 

where a is a positive parameter. 

While we were mainly motivated by trying to extend the classical probabilis- 
tic phenomena to the setting of free probability, it is worth mentioning that the 
theory of free extreme values is directly related to natural operations on random 
matrices (see the recent preprint [4]). The results of this paper can easily be 
translated in the context of [4]. 

The rest of the paper is organized as follows. Section [2] gives a brief intro- 
duction to free probability theory. Section [3] proves Theorem 1 using Theorem 
2. Section [4] details the definition of the convergence of free point process and 
proves Theorem 2. And Section [5] proves Theorem 3. 



2 Preliminaries 

2.1 Free Independence 

Definition 4 A W^* -probability space is a pair (A, E) , where A is a von Neu- 
mann algebra of bounded linear operators acting on elements of a complex sep- 
arable Hilbert space and E is a faithful normal trace that satisfies the condition 
E{I) — 1. Operators affiliated with algebra A are called non-commutative ran- 
dom variables, or simply random variables, and the functional E is called the 
expectation. 
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If P (dX) is the spectral resolution associated with a normal operator A, then 
we can define a measure fi{dX) = E {P (dX)) . It is easy to check that is a 
probability measure supported on the spectrum of A. We call this measure, /i, 
the spectral probability measure associated with operator A and expectation E. 

The most important concept in free probability theory is that of free inde- 
pendence of non-commutative random variables. Let a set of r.v. Ai, . . . , An be 
given. With each of them we can associate an algebra Ai, which is generated by 
Ai and A*; that is, it is the weak topology closure of all polynomials in variables 
Ai and A* . Let Ai denote an arbitrary element of algebra Ai . 

Definition 5 The algebras Ai, . . . ,An (and variables Ai, . . . , A^ that generate 
them) are said to be freely independent or free, if the following condition holds: 

E (A(l) . . . ^i(m)) = 0, 

provided that E (Ajj^-j) — and i{s + 1) 7^ i{s) for every s. 

For more information about non-commutative probability spaces and free 
operators we refer the reader to Sections 2.2 - 2.5 in the book [TB] by Voiculescu, 
Dykema and Nica. 

If X and Y are two free self-adjoint random variables with spectral proba- 
bilities measures fi and v respectively, then we denote the spectral probability 
measure oi X + Y a.s /j-Wv, and call it the free additive convolution of and v. 



2.2 Free Poisson random variables 

Let X be a self-adjoint operator that has the so-called free Poisson distribu- 
tion with parameter ("intensity") A. The continuous part of this distribution is 



supported on the interval 



1 - 



Va) ,(i + Vx 

\Jax-(1- X + x) 

l-KX 



and the density is 



In addition, if A < 1, then there is also an atom at zero with the probability 
weight 1 — A. We call such an operator X a (non-commutative) Poisson random 
variable with intensity A and size 1. 

The sum of two freely independent Poisson random variables of intensities Ai 
and A2 is again a Poisson random variable of intensity Ai -I- A2 (see, for example, 
a remark on pagelOS in [TO])- 

If we scale a non-commutative Poisson random variable by a, then we get a 
variable, which we call a scaled (non- commutative) Poisson random variable of 
intensity A and size a. 
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Non-commutative Poisson random variables arise when we convolve a large 
number, N, of Bernoulli distributions that put probability X/N on 1 and prob- 
ability 1 — X/N on 0. The following result is well-known, see (TU], [Il],or [TS]. 

Proposition 1 Suppose (n — 1,2,...) is a sequence of Bernoulli distribu- 
tions, such that fin ({!}) ~ X/n and fin ({0}) — 1 — fJ-n ({!}) • Define Vn as 
follows: 

I'ri = ffl • • • ffl Mn ■ 
^^^^^^^ 

n times 

Then Vn weakly converges to the free Poisson distribution with intensity X and 
size 1. 

2.3 Free Poisson random measure 
Definition 6 Let {0,B,h') be a measure space, and put 

Bo ^ {B e B : i^{B) < oo}. 

Let further {A, E) be a W* -probability space, and let A+ denote the cone of pos- 
itive operators in A. Then a free Poisson random measure ('fPrmj on z^) 
with values in {A, E) is a mapping M : Bq A+, with the following properties: 

(i) For any set B in Bq, M (B) is a free Poisson variable with parameter v (B) . 

(ii) If r € N, and Bi, Br € Bq are disjoint, then M (Bi) , AI (Br) are free. 
(Hi) Ifr e N, and Bi, Br G Bq are disjoint, then M (Uj^ii?j) = M (Bj) 

The existence of a free Poisson measure for arbitrary spaces (0,S, i^) and 
{A, E) was shown in [15] and a different proof was given in '2\ . 

Let / be a real-valued simple function in L^ (O, B, v) , i.e, suppose that it 
can be written as 

r 
i=l 

for a system of disjoint Bi E Bq. Then we define the integral of / with respect 
to a Poisson random measure M as follows: 

f f dM = J2 [B,] . 

It is possible to check that this definition is consistent. Moreover, as it is shown 
in [2], this concept can be extended to a larger class of functions: 

Proposition 2 Let f be a real-valued function in L^ {QjBjV) and suppose 
that Sn is a sequence of real valued simple B-measurable functions, satisfying 
the condition that there exists a positive v-integrable function h (6) , such that 
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\sn (0)1 l£ h (0) for all n and 0. Suppose also that lim„^oo Sn (0) — f (0) for all 9 
Then integrals Jq s„ dM are well-defined and converge in probability to a self- 
adjoint (possibly unbounded) operator I (/) affiliated with A. Furthermore, the 
limit I (/) is independent of the choice of approximating sequence s„ of simple 
functions. 

The resulting functional / (/) is defined for all real valued functions / in 
(Q,B,u) and is called the integral with respect to the free Poisson random 

measure M. It possesses all the usual properties of the integral: additivity, 

linear scaling, continuity, etc. 

3 Proof of Theorem 1 

As was noted in Introduction, only the equivalence of (i) and (iv) needs a proof. 
The equivalence of (i) and (iv) can be reduced to a problem about convergence 
of free point processes. Indeed, let /i„ (A) = ^ (a„A + 6„) . Then (i) is equivalent 
to the statement that (A) — > A (A) for all Borel sets Ac {x,oo) . 

Indeed, suppose that fi is in the domain of attaction of G {x) , and let F (x) 
denote the distribution function of the measure /i. Then 

F" (a„x + 6„)^G(x), 
For every x G {x, oo) , G (x) is positive, hence we can take logarithms and get 
n log F {unx + 6„) ^ log G {x) , 

which is equivalent to 

n{l -F{anX + bn)) ^-\ogG{x) = A((.x,oo)) . 

Consequently, 

n^n ((x, oo)) A ((x, oo)) , 

from which we conclude that n/i„ (A) A (A) for all Borcl sets A C (x, oo) . 

By reversing the steps of this argument we obtain the reverse implication: 
If (A) A (A) for all Borel sets A C (x, oo) , then fi is in the domain of 
attaction of G (x), and (i) holds. 

Therefore the equivalence of (i) and (iv) follows from Theorem [2] if we take 
{Xi — On) /bn as the triangular array Xi^„- 

4 Proof of Theorem 2 

4.1 Weak Convergence 

In this section, we define precisely the mode of convergence of free point mea- 
sures that we use. It corresponds to the weak convergence of point processes in 
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the classical case. 

Let D he a Borel subset of M and let (D) denote the space of bounded, 
Borel measurable functions that have compact support on D. 

Definition 7 We say that a free point process Mn converges weakly on D to o 
free Poisson random measure M, which is defined on {D,B, A) and takes values 
in A, if for every function f G J-'^ (D) the following convergence holds: 



Sometimes we also need to speak about convergence with respect to a class 
of functions, which is different from JFJ^ (D). 

Definition 8 We say that a free point process Mn converges weakly with re- 
spect to a class of functions J- to a free Poisson random measure M, if for every 
function f € J- the following convergence holds: 



We will prove Theorem [2] by considering initially the convergence of free 
point processes M„ with respect to the class of simple functions (i.e., finite 
sums of indicator functions), and then approximating functions from a more 
general class by simple functions. 

4.2 Convergence with respect to simple functions 

Let S (D) be the class of simple functions on C M, i.e., the class of finite sums 
of indicator functions of Borel sets belonging to D. 

Proposition 3 Let Xi n be a triangular array of free, self-adjoint random vari- 
ables and let the spectral probability measure of Xi^n be fJ-n- Let X be a Radon 
measure on D QR. If 



for each Borel set A d D, then the free point process M„ associated with the 
array „ converges weakly with respect to S (D) to a free Poisson random 
measure M with the intensity measure A. 

Before proving this proposition, we derive some auxiliary results. 

Lemma 1 Suppose Xi^„ is an array of free and identically distributed random 
variables with the spectral measure pin- Let npLn {A) — > X{A) < oo as n oo. 
Let Zi,n = {Xi,n) ■ Then as n — > oo, the sum Sn = X]"=i ^i-n converges in 
distribution to a free Poisson random variable with intensity A (A) . 





nun (A) A (A) 
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Proof: Note that Zi.n are projections with expectation (A) and they 
are free. Therefore, X^iLi -^i," ^^"^ sum of free projections and we can use 
Proposition [T] to infer the claim of the lemma. QED. 

As the next step to the proof of Proposition [3] we need to check that if Borel 
sets Afc are disjoint, then the sums — X^iLi l'4fc (^i-n) ^^'S asymptotically 
free with respect to growing n. 

Recall the definition of the asymptotic freeness: Let [Ai, Ei) be a sequence of 
non-commutative probability spaces and let Xi and Yi be two random variables 
in Ai- Let also x and y be two free operators in a non-commutative probability 
space {A, E). 

Definition 9 The sequences Xi and Yi are called asymptotically free ij the 
sequence of pairs {Xi,Yi) converges in distribution to the pair {x,y) . That is, 
for every e > and every sequence of k-tuples (ni,...,nfc) with non-negative 
integers nj, there exists such io that for i > ia, the following inequality holds: 

\E, - £;(a;"iy"^..a;"'=-iy"'')| < e. 

At the cost of more complicated notation, this definiton can be generalized 
to the case of more than two variables. 

Lemma 2 Let P^^, (where n = 1,2, i ~ 1, n, and fc = 1, be pro- 
jections of dimension A^'^^' /n. Assume that for each n, algebras Ai generated by 
sets \Pin\ '^'^^ free. Also assume that for each n and i, the projections 

^i'n ^'''^ orthogonal to each other, i.e., Pi^^Pin for every pair k ^ k' . Let 
sli^ = X]r=i Pi'n ■ Then as n ^ oo, the sequences S*!*^^ converge in distribution 
to freely independent variables 5^'"'^ that have free Poisson distributions with pa- 
rameters A''^\ respectively. In particular, the sequences are asymptotically 
free with respect to growth in n. 

Proof: The fact that each of the sequences Sn converge in distribution 
to a variable S'^'^^ that has a free Poisson distribution is clear from Proposition 
[TJ The essential part is to prove that asymptotic freeness holds. This claim 
is a direct consequence of Speicher's multidimensional limit theorem (see, for 
example. Theorem 13.1 in the book [12] by Nica and Speicher). Indeed, we 
need to prove that all mixed free cumulants of the limit are zero. By Speicher's 
theorem, this is equivalent to the statement that the following limits are zero: 

iimni?(<^'p(tr)...p(t:')=o. 

Here ki, . . . ,ks is an arbitrary s-tuple with the property that it has a pair of 
distinct coordinates, i.e., ki ^ kj. However, the fact that these limits are zero 
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is clear from the assumption that the projections J are orthogonal to each 
other. QED. 

Now we can proceed to the proof of Proposition [3l 

Proof: Let / = X]fc=i Cfcl^fc (x) , where Ak are disjoint Borel sets. Using the 
assumption that n/z„ (Ak) A (Ak) and Lemma [1] we can find a free Poisson 
random measure M such that 

n „ 

V (X,.„) 4 Af (Afe) = / (a;)M(dx) 

Jr 

as n —f oo. Indeed, it is enough to take a Poisson random measure M with the 
intensity measure A. 

In addition, by Lemma[2l sums Sk = ^Af. (Xi^n) become asymptotically 

free for different fc as n grows. Since M (Ak) are free by the definition of the 
free Poisson measure, this implies that 

r n r ^ r 

Y.CkY. I A, 4 J2 CkM {Ak)^J2 / 1-4, {x) M {dx) . 

k=l i=l k=l fc=l "^^ 

as n — !■ cx). Therefore, 

n „ 

Y,fiX^,n)^ / f{x)M{dx), 
^=l 

where we used the additivity property of the integral with respect to a free 
Poisson random measure (see [5], Remark 4.2(b)). QED. 



4.3 Convergence with respect to bounded, Borel measur- 
able functions with compact support 

The goal of this section is to prove our main Theorem [2l 

Consider a bounded, Borel measurable, compactly supported function / : 
Z? ^ R, such that < / < 1. (A more general case of a function /, which satisfies 
Ci < / < 6*2, can be treated similarly.) For positive integers N = 1,2,..., and 
fc = 1, . . . , iV, define the set 

4"^ = {-eBupp(/):^</(.)<A|. 

The sets A^j!^"^ are disjoint, measurable, and have finite A-measure. Their union 
is D. 

We define lower and upper approximations to the function / as follows: 

k=l 
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and 

A k 

k=l 

We note that: 

(i) (x) < (x) ; 

(ii) l'^ (x) is an increasing sequence of functions; 

(iii) (x) is a decreasing sequence of functions, and 
iv) liniAr_^oo (x) — {x) = uniformly in x. 

The functions {x) and u'^ {x) are simple: (x) — X^t^i (^) 

and (x) = 4^^1a<"' (^) ' ^o*^ ^1^° ^^Pt (4^^ " " 

converges to zero as A'^ ^ oo. 

Let us drop for convenience the superscript N when we consider it as fixed, 
and simply write / (x) = J2i=i ^k^Ak {x) and u {x) = X]i=i dklA^ {x) , where Ak 
are disjoint Borel-measurable sets. By Proposition [3l as n — > oo, 

n N 

1=1 fc=i 

where Mj, are freely independent Poisson random variables with intensities = 
A (Ak). Let Fi (x) denote the distribution function of J2k=i c/c-^fe- 
Similarly, 

n N 

t=i fe=i 

and we denote the distribution function of X^feLi dkMk as Fu (x). 

Let Ff^n denote the distribution function of f i^i.n) and let Ff be one 
of the limit points of this sequence of distribution functions. 

Proposition 4 Ff is a distribution function and F^ {x) < Ff (x) < Fi (x) for 

every x. 

Proof: We will infer this from Lemma [3] below and its Corollary. This 
lemma is a particular case of Weyl's eigenvalue inequalities for operators in a 
von Neumann algebra of type Hi . If Fj^ (x) is the spectral distribution func- 
tion of a self-adjoint operator A, then we define the eigenvalue function 6*^ (t) = 
inf {x : Fa{x) > 1 — t} . The function 9a {t) is non-increasing and right-continuous. 
Intuitively, it can be thought of as a "sequence of eigenvalues" of A, indexed in 
decreasing order by parameter t. 

Let us use notation 9a {t — 0) to denote lim^j^o 9a {t — e) . Then the following 
generalization of Weyl inequalities holds: 
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Lemma 3 If A and B are two hounded self-adjoint operators from a W*- 
probability space A and if B is non-negative definite, then 

OA{t) < eA+B{t)<eA{t) + \\B\\, and 

0A{t-O) < 9A+Bit~0)<eA{t-0) + \\B\\. 

Corollary 1 If B > 0, then (J-a+b ^ /^A, that is, Fa+b (x) < Fa (x) for each 

X. 

Proof of Lemma [3j These results easily follow from an inequality in [5] 
which states that if (a — e,a) C [0, 1], (6 — e, b) C [0, 1] , and a + b < 1, then 

ra-\-b pa pb 

/ Oa+b (t) dt< 9 A (t) dt+ 9b (t) dt. (2) 

'J a-\-b—e J a—e J b—e 

QED. 

By Corollary [U for each n the distribution Ff^n is between the distribution 
functions of Y^^=i " {^i,n) and ' {^i,n) ■ As n grows, these two sequences of 
distribution functions approach Fu {x) and Fi {x) , respectively. Therefore, every 
limit point of -F/,n is between F^ and Fi. The claim that Ff is a distribution 
function follows from the fact that both F^ and Fi are distribution functions. 
QED. 

Now we want to show that i^i^' {x) approaches f/^"* (x) as N grows. 
Recall that the Levy distance between two distribution functions is defined 
as follows: 

dL {Fa, Fb) = supinf {s > Q : Fb {x ~ s) - s < Fa{x) < Fb {x + s) + .s } . 

We can interpret this distance geometrically. Let T a be the graph of function 
Fa, and at the points of discontinuity let us connect the left and right limits by a 
(vertical) straight line interval. Call the resulting curve T a- Similarly define V b- 
Let d be the maximum distance between and in the direction from the 
south-east to the north-west, i.e., in the direction which is obtained by rotating 
the vertical direction by 7r/4 counter-clockwise. Then dL {Fa, Fb) = d/^/2. 

Proposition 5 Let K be the sum of intensities of freely independent Pois- 
son random variables and let Fi {x) and F^ {x) be distribution functions 
ofJ2k=i^kMk andJ2k=idkMk Then 

dL{Fi,Fu)<(2K + 3VK+l) sup {dk - Ck) ■ 

^ ' l<k<N 

Remark: In the proof of Theorem [21 the finiteness of K will be ensured by 
the assumptions that measure A is Radon and that / has a compact support. 
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For the proof of this proposition we need two lemmas. Lemma H] provides 
a bound on the norm of the sum of scaled Poisson random variables in terms 
of the sizes of these variables, and Lemma [S] relates the Levy distance between 
two random variables to the norm of their difference. 

Lemma 4 Let Mi, (i — 1, ...,rj be freely independent Poisson random variables, 
which have intensities A.;, and let bi be non-negative real numbers. Assume that 
X]i=i — ^ '^'^'^ ^ ~ sup]^<;,;<^ bi. Then 



< b (2K + sVk + 1 



Proof: Let Xi be free self-adjoint random variables that have zero mean. 
Then by an inequality from [17j : 



< max \\Xi\ 

l<i<r 



J2Var{X, 



If Yi are free self-adjoint random variables with non-zero mean, and Xi = Yi 
E {Yi) , then the previous inequality implies that 



< 



< 



r 



i=l 



max d (Yi) 

l<i<r 



Y.Var{Y{), (3) 
\ '=1 



where d (Yi) is the diameter of the support of Yi. 

We will apply this inequality to Yi — biMi and estimate each of the three 
terms on the right-hand side of ([3]) in turn: 

1) Since E {Mi) = Xi, and J2>'i<K, therefore XILi (^^0 < bK. 

2) The diameter of the support of biMi is less or equal to bi (l + \/Xi)^ < 
b(l + 2VK + K] 



3) Since Var (Mi) = Xi, therefore VELi (biM,) < b^K. 
In sum, biM^W < b (2K + 3\/X + l) . QED. 



Lemma 5 Let A and B be two bounded self-adjoint operators from a W* 
probability space A and assume that B — A>{). Then 



dL {Fa.Fb) < \\B-A\\ 
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Proof: Let Fa and Fb be distribution functions, and 9a and Ob be the 
corresponding eigenvalue functions. Then we claim that 

dL{FA,FB)< sup \eAit)-0B{t)\. (4) 

0<t<l 

Indeed, let the graphs of functions 6 a and 9b he denoted as Aa and Ab, respec- 
tively. Connecting the left and right limits at the points of discontinuity gives 
us the curves and A^- It is easy to see that these curves can be obtained 
from curves Ta and Tb (i.e., the graphs of Fa (x) and Fb (x) with connected 
limits at the points of discontinuity) by rotating them around the point (0, 1) 
counter-clockwise by the angle tt/2 and then shifting the result of the rotation 
by vector (0, —1) . It follows that the distance d, which was used in the defini- 
tion of the Levy distance can also be defined as the maximum distance between 
Aa and Ab in the direction from the south-west to the north-east, i.e., in the 
direction which is obtained by rotating the vertical direction by 7r/4 clockwise. 
Since 9a (t) and 9b {t) are non- increasing functions, therefore 

d<V2 sup \9A{t)-9B{t)\. 

0<t<l 

This implies {Fa,Fb) < supo<t<i \9a (t) - 9b (i)| • 

Inequality ([U and Lemma [3] imply the statement of the lemma. QED. 
Now we can prove Proposition (5] 

Proof of Proposition [5j Let X = J2k=i i^k — Ck) Mk- By Lemma HI 
\\X\\ < b (2K + 3\/K + , where b — supj^^j,^^ (dfc — Cfe) and K is the sum of 

the intensities of Mk- By LemmaO this implies that cIl (Fi, F^) < b (2K + 3\/K + 1 
QED. 

Using Proposition m we can proceed to the proof of Theorem [2l By Propo- 
sition [Sj we know that if is fixed and n — > oo, then 

n N 

i=l k=l 

and 

n N 

4=1 k=l 

where M is a free Poisson random measure with intensity A (dx) . Let the dis- 
tributions of the right-hand sides be denoted as Fin and F^^n . 

By Corollary [1] (p. [T5)) . Fjn is a decreasing sequence and i^„jv is an increasing 
sequence of distribution functions. In addition, Fin (x) > F^n (x) for every N 
and X. Since the sum of intensities of variables M ^^[.^''^ is less than A {D) < 
oo by assumption, therefore Proposition [5] is applicable and we can conclude 
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that the Levy distance between Fin and F.^^n converges to zero as TV cx3. 
Consequently, these two distributions (weakly) converge to a limit distribution 
function as N ^ oo. 

Moreover, by the definition of the integral with respect to a free Poisson 
random measure, this limit equals the distribution function of / / {x) M [dx) . 

In addition, by Proposition |4] every limit point of the sequence of Ff^n is 
between Fin and F^n for every TV, and therefore the sequence of Ff^n also 
converges to the distribution function of / / {x) M (dx) . QED. 

This completes the proof of Theorem [2] 

5 Proof of Theorem 3 

Recall that we defined the distribution of a free order statistic in the following 
way. Let Xi, ...,X„ be freely independent self-adjoint random variables and let 
Xi have the spectral distribution functions Fj. Let 

n 
1=1 

Definition 2. For every real k > 0, we say that Fn.k (t) =: E [l[o.fc] {Yn (t))] 
is the distribution function of k-th order statistic of the sequence Xi, . . . ,X„, 
and that it is the k-th order free extremal convolution of distribution functions 

It is straightforward to check that in the case of commutative random vari- 
ables, this definition gives the distribution function of the usual {\k\ -\- l)-order 
statistic. 

In the non-commutative case, we need to check that this is a consistent 
definition, and that Fn^k (t) is indeed a probability distribution function for 
each A: > 0. 

It is easy to see that Fn^k (t) is non-decreasing in t. Indeed, let t' > t. Then 
for each i, l(f/.oo) i^i) < l(t,oo) i^i): and therefore, Y {t') < Y (t) . It follows 
that l[o,fc] (Y it')) > l[o,fe] (Y (t)) , and therefore F„,fc (f) > F^^k (t) . 

This function is also right-continuous in t. Consider a sequence tm I t. First, 
note that l(i,„^oo) i^i) ~^ l(i,oo) i^i)- Second, since operators 1(4, oo) i^i) are 
freely independent for diffferent i, this implies that Y (tm) Y (t) as tm [ t. 
Indeed, the operators Y (tm) and Y (t) are uniformly bounded (||i^(tm)|| < n 
and < n), and the moments of the distribution of Y (tm) converge to 

the corresponding moments of the distribution of Y (t) . 

Third, let the spectral probability distribution functions of Y (tm) and Y (t) 
be denoted as Gm (x) and G (x) , respectively. Then E [l[o,fc] (Y (<m))] — Gm (k) 
and E [l[o,fc] (Y (t))] = G (k) . Since G™ (k) = F„,k (i™) , and G (k) = F„,k (t) , 
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therefore we aim to prove that Gm (k) G {k) as m cx3 for all k. The con- 
vergence Y (tm) — * Y (i) means the convergence of the moments of the spectral 
probability measures of operators Y (tm) and Y (t) , which implies weak conver- 
gence of these measures because the measures have uniformly bounded support. 
This implies that Gm (k) — > G (fc) as m ^ cx), for all points k at which the prob- 
ability distribution function G (k) is continuous. We will prove that, moreover, 
even if G {x) has a jump at a; = fc, then the sequence (fc) still converges to 
G (fc) . At this point of the argument, it is essential that tm converges to t from 
above and therefore Gm {k) > G (k) . 

Indeed, by seeking a contradiction, suppose that Gm (fc) does not converge 
to G (k) . Then, take e such that (i) Gm {k) — G{k) > e for all m, and take k' > k 
such that (ii) k' is a point of continuity of G (x) , and (iii) G (fc') — G (k) < e/2. 
Such k' exists because G (x) is a spectral probability distribution function and 
therefore it is right-continuous. Since Gm (k) is increasing, we conclude from 
(i), (ii), and (iii) that Gm{k') — G [k') > e/2 for all m. But this means that 
Gm (x) does not converge to G (x) at a point of continuity of G (x) , namely, at 
k' . This is a contradiction, and we conclude that Gm (k) converges to G (fc) for 
all k. This means that (t) is right-continuous in t. 

Finally, as i ^ oo, l(^t,oo) 0. Therefore Y {t) 4 0, and 1[om] {Y {*)) ^ 

/. Hence F„ ^ (t) ^ 1 as t ^ oo, and we conclude that Fn^k (t) is a valid 
distribution function. 

Consider now the special case when fc = 0. In this case Fnfl (t) is the dimen- 
sion of the nill-space of Yn (t) , which equals to the dimension of the intersection 
of the nill-spaces of l(t,oo) i^i) ■ It is easy to see that this coincides with the 
definition of the free extremal convolution of the distribution F, which was in- 
troduced in [3]. 

Now let us investigate the question of the limiting behavior of the distribu- 
tions F„ k (t) when n — *■ oo. The limits are described in Theorem [31 
Proof of Theorem [3} For each n we re-define: 



where M„ is the free point process associated with the triangular array (Xi — hn) /a. 



which is a free Poisson random variable with the intensity A (t) = — log G {t) . 
Then, in order to calculate the limit of Fn^k (t) for n — > oo, we only need to 




n 



i=l 



The bracket (M„, l(t,oo)) 



converges in distribution to a random variable Gf 



calculate i?l[o,fe] (Gt) , that is, the distribution function of Ct at k. Let us denote 
the distribution function of Ct as Ht (x) , 
For fc < 0, we have Ht (fc) = 0. For fc = 0, 
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For fc > 0, 



Ht (k) = { 



where 



Ht (0) 



Ht (0) , 

1 



2 Pt (0 if ^ e 



if A: < 1 - 



if fc > 



, 1 



^4C-(l-A(t) + 0^ 

2< 



Then, we need to compute i^jj.) (t) , which is Ht (fc) considered a function of t 
for a fixed k. Let (cc) denote the solution of the equation X{t) = x. (That is, 
if [x) is the functional inversion of the limit distribution function G {t) , 

then A-i [x] = G'(-i) (e"^) .) 

Then, for k = 0: 



0, ifi<A-i(l), 
l-A(t), ift>A-i(l). 



For fc e (0,1): 



Fik) it) = { 



0, a t < x-^ y(i + Vk^ 

J'l^.^yPtiOd^, ifie A-1 (^(l + Vfc)'),A-i(l) 

l-^W + /(i_^)=Pt(Orf^, if (tGA-Ml),A-^((l-Vfc)') 

l-A(t), ifi> A-i . 

For k > 1, we have: 



0, 



^'(fc) it) = { 



if i e 



if t < A-1 (^(l + VkY^ , 

i((i + yfc)'),A-Mi) 



1 - ^ (0 + (0 if (t e A-i (1) , A-1 (^(l - A/fc' 



if i > A" 
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Fik) it) 



Combinining these cases, we obtain the foUowing equation: 

0, iit<X-^ 

l-A(i)l[o,i)(A:), ifi> A^i (^(l-Vfc) 

QED. 



l + Vk 
(^(l + x/fc)'),A-i(l) 



Example 4 Distributions from the domain of attraction of Type II extremal 
value law 

Consider the case of convergence to the Type II extremal vahie law, when 
the constants a„ and 6„ are chosen in such a way, that the limit law is G (x) = 
exp {—x~'^) for X > O.Then we can conclude that the limit distribution of the k 
order statistic is given as follows: 



if < < I 1 + Vfc 
if t e 



1 + ,1 



i-t-'' + I^,_t^.f2yPti0dC, ate i^i,i^{i-Vk) 
1 - t-n[o,i) (k) , ift>((l-Vk 



2\ -l/"^ ■ 



where 



4e-(i-i- + e)' 



2< 



We illustrate this result for some particular values of i' and k. 
Consider fc = 0. Then 



0, 

1 - r 



if i < 1, 

if i > 1. 



This is the Type 2 ("Pareto") limit distribution in Definition 6.8 of [3]. 

The distributions of k order statistics for different values of k are illustrated 
in Figure 1. 



[Put Figure 1 here.] 
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It is interesting to note that if A; > 1, then for all sufficiently large t, F(^i^) (t) = 
1. This can be interpreted as saying that the scaled k order statistic is guaranteed 
to be less than for a sufhiciently large Iq. In another interpretation, this result 
means that for our choice of scaling parameters a„ and 6„ and for every fc > 1, 
if t is sufficiently large, then 



for all large n. 

A similar situation occurs in the classical case if the initial distribution (i.e. 
the distribution of Xi) is bounded from above. In this case the limit distribu- 
tion is also bounded from above. In contrast, in the free probability case this 
situation occurs even if the initial distribution is unbounded from above. Our 
previous example shows that this situation occurs even if the initial distribution 
has heavy tails. 
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